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Abstract 

We consider the first-crossing-time problem through a constant boundary for a 
Wiener process perturbed by random jumps driven by a counting process. On the base 
of a sample-path analysis of the jump-diffusion process we obtain explicit lower bounds 
for the first-crossing-time density and for the first-crossing-time distribution function. 
In the case of the distribution function, the bound is improved by use of processes 
comparison based on the usual stochastic order. The special case of constant jumps 
driven by a Poisson process is thoroughly discussed. 
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1 Introduction 

In a variety of applied contexts a relevant role is played by jump-diffusion processes, i.e. by 
diffusion processes to which jumps occurring at random times are superimposed. Indeed, 
such processes are for example invoked for the description of stochastic neuronal activity 
(see Giraudo and Sacerdote [13^ and Giraudo et al. [T5|), of complex queueing systems (see 
Perry and Stadje tlZ|)j of random assets in mathematical finance (see Ball and Roma [4]), 
of surplus of insurance companies in ruin theory (see Gerber and Landry [llj). of acto- 
myosin interaction in biomathematics (see Buonocore et al. [6]). Despite the relevance of 
the first-crossing-time (FCT) problem for jump-diffusion processes in such contexts, only few 
analytical results are available on their probability density function (pdf), even in the case 
of very simple boundaries. The available results are mainly focused on equations involving 
the FCT moments (see, for instance, Abundo [2], Giraudo and Sacerdote [T^], and Tuckwell 
|21j). which however seem to be hardly manageable for practical purposes. Other results 
concerning bounds, obtained by use of Laplace transform, are limited to upper and lower 
bounds for the ruin probability in a jump-diffusion process involved in a risk model perturbed 
by Brownian motion (see Yin and Chiu [22]), and for the mean and the variance of the hitting 
time in certain jump-diffusion processes (see Schal [T^). Analytical results on FCT pdf's are 
also very rare (see Kou and Wang J^, where explicit solutions of the Laplace transform of 
the distributions of the first-crossing times are disclosed for a Brownian motion perturbed by 
double exponentially distributed jumps). Hence, as a viable alternative, efficient algorithms 
have been devised in order to evaluate FCT densities (cf. the recent contributions by Atiya 
and Metwally [3], and by Di Crescenzo et al. [7 ). 

In order to evaluate the performance of simulation algorithms, one needs to come up with 
some sample cases in which FCT densities and distribution functions obtained by simulation 
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can be compared with the corresponding bounds analytically determined. Hence, in the 
present paper lower bounds for FCT densities and distribution functions will be determined 
for jump-diffusion processes based on the Wiener process in the presence of a constant 
boundary. Our results are based on a sample-path analysis of the jump-diffusion process 
and on some specific features of the underlying Wiener process, such as the space and time 
homogeneity and the availability of a closed form of the FCT density through a constant 
boundary. It must be pointed out that our approach is at all different from that of Bischoff 
and Hashorva [5], where the Cameron-Martin-Girsanov formula is exploited to obtain a 
lower bound for the boundary crossing probability of Brownian bridge with trend. 

In Section 2 we formally describe the jump- Wiener model, that consists of the superpo- 
sition of a Wiener process and of a jump process with generally distributed jumps occurring 
at the occurrences of a counting process. The FCT problem for such process through a 
constant boundary is then addressed and a lower bound for the FCT pdf is then determined 
in Section 3. Use of such a bound is then made for the special case when upward and down- 
ward jumps have constant amplitudes and occur according to a Poisson process. Section 4 
presents a lower bound for the FCT cumulative distribution function (cdf). The special 
case when only upward constant jumps are allowed, and are separed by random times hav- 
ing exponential distribution is thoroughly investigated. In this case the bound is improved 
by making use of a technique based on the comparison of FCT's by the "usual stochastic 
order" . Finally, in Section 5 some remarks on the computational aspects are given. 

2 FCT problem for a jump-diffusion process 

Let {X{t)}t>o be a jump-diffusion process defined by 



where {W{t)}t>Q and {Y{t)}t>o are independent stochastic processes and 

(i) {W{t)}t>o is a Wiener process with drift /i e K, and variance E (0, +oo) per unit 

time, starting at W{0) = xq; 

(ii) {y{t)}t>Q is a jump process such that Y{0) = and Y{t) — J^i'Ji --^i^ ^ > 0, with 

Y{t) = when N{t) = and Ji, J2, . . . real- valued i.i.d. r.v.'s such that V{Ji = 0} < I 
for all i. The cdf of J,j will be denoted by Fj{x) and the survival function by Fj{x) = 
1 — Fj{x). By {N{t)}t>o we denote a counting process independent of { Ji, J2, . . .} 
and characterized by i.i.d. absolutely continuous positive renewals i?2, • ■ • having 
pdf /rIx), cdf Ffi{x) and survival function Fj^{x) = 1 — Fj^{x). 



(1) 



X{t) = W{t)+Y{t) 



Let 



fix, t\xo) = ^P {X{t) < X I X{0) - xo} , 
be the conditional pdf of {X{t)}. Since 



t > 0, 




where denotes the cdf of J2'i=i Ji 
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is the pdf of W(t), from the assumptions on {X{t)} we have: 

(2) /(a:,t|a:o) -^/'[A^(i) = fc] / fw{x-u,t\xo)dF^''\u). 

Denoting by (j)iq{u,t) the probabihty generating function of N{t) and by ipj{s) the moment 
generating function (m.g.f.) of J,;, the m.g.f. of {X{t)} is given by 

(3) £; [e"^(*) I X(0) = xo] =0w[V'j(s),i] exp|(a;o + /xi)s + ^s2 
and 

(4) E[X{t) I X{0) = xo] = X0 + + E{.h)E[N{t)l 

(5) Ya.v[X{t) I X{Q) = xo] = a^t + Var(Ji) E[N{t)] + £;2(Ji) Var[iV(t)]. 

Assuming xq < S, hereafter we shall discuss the FCT problem through a constant 
boundary S for the jump-diffusion process defined in |T]). We denote by 

(6) Tx = inf{t > : X{t) > S}, P{X{0) = xo} - 1, 
the FCT of {X{t)} through S from below, and by 

(7) gxiS,t\xo)^-^^P{Tx<t\X{0) = xo}, <>0 
its pdf. The corresponding cdf will be denoted by 

(8) Gx{S,t\xo)=p{Tx <t\X{0) = xo), t>0. 

To determine the ultimate FCT probability and the FCT moments one is usually asked to 
solve appropriate integro-differential equations (see Abundo [2] and Tuckwell HJ). The de- 
termination of closed-form expressions for the FCT pdf and cdf is instead a harder problem, 
because so far no analytical method appear to be available thus for. Hence, in the following 
sections we shall confine our investigation to determining useful lower bounds for density 
© and cdf (0). 

Let Tw denote the FCT from below of Wiener process {VF(t)} from xq < S to the 
constant boundary S. Some well-known results on Tw that will be used later are recalled 
hereafter: 

• FCT cdf of {W{t)} through S: 

(9) Gw{S, t\xo) = P{Tw < 1 1 W{0) - xo} 

/ S~xo-nt\ , fc, S~xq\ f S-XQ+^it 
= * r=T7 + exp 2^ — $ 



where ^ 
(10) $(z) = / e-^'/^dx, z e 



/27r J^c 

• FCT pdf of {W{t)} through S: 

(11) g^iS,t\xo)^i-P{Tw<t\W{0)^xo}^^=^e^p-^ {S-x.-^tf 



dt ' V2TraH^ { 2aH 
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iS'-avoiding transition pdf of {W^(t)}: 



(12) awixM^o) = Q^V{W{t)<x,Tw>t\W{Q)^xa} 



1 { [x-xf)- fity 

cxp ■ 



^ 2{S-x)iS-xn) 
^-^"P^ 



V2t^ I 2<jH 

3 Lower bound for FCT pdf 

In order to obtain the preannounced lower bound for pdf ([7]) let us note that, for alH > 0, 
event {Tx S {t,t + dt)} can be decomposed into the following three mutually exclusive 
events: 

(i) the first jump occurs after t and the first crossing through S occurs in (t, t + dt) due to 
the diffusive component of {X{t)}; 

(a) in (0,t) the diffusive component of {X{t)} does not cross the boundary, the first jump 
occurs in (t, t + dt) and it causes the first crossing; 

(in) the first jump occurs at time 6 £ (0,i) and it does not cause the first crossing, the 
diffusive component of {X{t)} having not crossed the boundary in (0,9); the first crossing 
finally occurs in {t, t + dt). 

Hence, for alH > the following equation holds: 

(13) gx{S,t\xo) = Ffiit) gwiS,t\xo) + fnit) / aw{x,t \ xq) Fj{S - x) dx 



t nS-x 

dFR{9) / aw{x,9\xQ) / gxiS,t - 9 \ x + u) dFj{u) dx. 

J —oo J —oo 

The formal proof of Eq. (fT3| has been given by Giraudo and Sacerdote [M] in the special case 
of jumps of constant size at the occurrence of a Poisson process, aiming to find out conditions 
under which the FCT density becomes multimodal (see also Sacerdote and Sirovich [T5] on 
this topic). As an immediate consequence of (|T3| a lower bound for the FCT pdf of {X{t)} 
can be obtained. Indeed, for alH > the following inequality holds: 

(14) gx{S,t\xo)>FR{t)gw{S,t\xo) + fRit)f awix,t\xo)Fj{S-x)dx 

J —OO 

+ FR{t-9)dFR{9) aw{x,9\xQ) gw{S,t ^ 9\x + u)dFj{u)dx. 

Jo J —oo J —oo 

Note that a tighter lower bound can be obtained by repeated substitutions of gx in the 
last term of the right-hand-side of l|13p . However, this would include many terms involving 
progressively high-order integrals, unsuitable for computational purposes. 



3.1 A special case 

Let us now study a special case of model ([T]). First of all, we assume that the jumps 
are separated by i.i.d. exponential random times Ri, i.e. restrict our attention to the case 
when {N{t)} is a Poisson process whose parameter will be denoted by A. In this case pdf 
f{x,t\ xq) is solution of the following integro-differential equation (see Buonocore et al. [B] 
or Di Crescenzo et al. [9]): 
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which clearly reduces to the Fokker-Plank equation of a Wiener process if A = 0. Moreover, 
let us assume that the jumps have positive or negative constant amplitude. In other words, 
the random jumps Ji are distributed as 



(15) 



a w.p. r] 
~b w.p. 1 — ?7, 



with 0<77<l,a>0 and 6 > 0. Hence, Y{t) can be expressed as 

Y{t) = aNi{t) ~bN2{t), t>0 

where {Ni{t)}t>o and {N2{t)}t>o are independent Poisson processes with rate 77A and (1 — 
r])X, respectively. From for x S R and t > we obtain 



(16) /(x,t|xo) 



-At 



00 00 



^ ^ (rjXty ((1 - 77)A^)fe ^^^^^ 



fc=0 j=0 



J! 



fc! 



(x + aj — bk — Xq — /ii)^ 
2^t 



Since under the present assumptions ipj{s) = 776°" + (1 — 77) e and (j)N{u,t) — e ^^'^'^ "\ 
from ^ the m.g.f. of {X{t)} for s € R becomes 



E 



X{Q) ^xo ^ exp <j -At [1 - r; e"^ - {I - v) + i^o + ^^t)s + 



As already pointed out in Di Cresccnzo et al. [71 , from ^ and ([5]) we have 

E[X{t) I X{Q) = xo] = xo + ^lt+ [yya - (1 - 7^)b]\t, 
Ya.v[X{t) I X{Q) = xo] = a^t + [r^a^ + (1 - ri)b^]\t. 

We stress that the lower bound given in (|14p can now be explicitly evaluated. Indeed, 
under the present assumptions for all < > we have: 



gx(S,t\xQ) > gi{S,t\xo), 



with 



(17) geiS,t\xo) -.^ e < gwiS,t\xQ) + t]X awix,t\xo)dx 

[ Js-a 

+ rjX d9 aw{x,6\xo) gw{S,t — 6\x + a)dx 



+ (1 — ?7)A d9 aw{x,0 \xo) gwiS,t ~ 9\ X — b) dx > , 



where gwiS, t\xo) is given in (fTTj) and 



aw{x, t \ xo) dx = ^ 



S-a 



S — xq — 



S — a — xq — lit 



-exp \^{S ~ xa) 



* - 



S ~ Xq + fit 

aVt 



S + a — Xq + fit 
aVt 



with $(2) defined in (fTO|) . The right-hand-side of (I17|) can be numerically evaluated. Indeed, 
making use of identity 



poo 

/ (z-5)e-("^ 
^0 



2a 



2a 



e4a 



1 - $ 



a > 0, 
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Figure 1: Estimated FCT density ^x{S, t\xo) for the Wiener process with constant Poisson- 
paced jumps for a ~ b ~ 3.75, and fower bound ge{S,t \ xq), with xq — 0, S — 5, X = 0.2, 
ry = 1, CT^ = 0.2 and (a) /i = 1.0, (b) /i = 1.5, (c) ^ ^ 2.0, (d) = 2.5. 



following from Eq. 7.4.2 of Abramowitz and Stegun Ij, one obtains: 

5 — m 



awix,9 \ xo) gwiS,t ~ 9 \ x + u) dx = 



1 



■ exp 
- exp 



(5 — u — xo — /it) 
2(5" — Xi^)u 



A+{9,u) 



where m = maxjit, 0} and 

[(to - u)9 ±{S-XQ±m){t- 9)f 



A±{9,u) — exp 



27r(t - t 
aH9 



2a^et{t - 9) 
(S* - xo ± u) < 1 - $ 



(to - u)9 ± (S* - a;o ± m){t - 9) 
^a^9t{t - 9) 



Figures 1 and 2 show certain multimodal estimates of gx{S, t \ xq) obtained by means of 
the simulation procedure described in Di Crescenzo et al. [J , together with the lower bound 
given in (fT7|) . The latter appears to be very close to the estimated pdf in various cases; the 
goodness of the approximation is discussed in Section 5. 



4 Lower bound for FCT cdf 

In this section we obtain a lower bound for the FCT distribution function ([8]). Similarly to 
the case of the FCT density we note that event {Tx < t}, t > 0, can be decomposed into 
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Figure 2: Estimated FCT density 'gx{S, 1 1 xq) for the Wiener process with constant Poisson- 
paced jumps, and lower bound gi{S,t\xo), for jj, = 2.0 and (a) -q = 0.8, (b) 77 = 0.5, (c) 
T] = 0.2, (d) rj = 0. All other parameters are chosen as in Figure 1. 



three mutually exclusive events: 

(i) the first jump occurs after t and the first crossing through S occurs in (0, t] due to the 
diffusive component of {X{t)}; 

(a) the first jump occurs at time 9 e (0,t) and it causes the first crossing through S\ 
(in) the first jump occurs at time 9 e (0, i) and it does not cause the first crossing, the 
diffusive component of having not crossed the boundary in (0,6'); the first crossing 

finally occurs in {9,i\. 

The following equation thus holds for all t > Q: 

(18) Gx{S,t\xQ) =FR{t)Gw{S,t\xo) + f dFai^) [ awix,^\ xo)Fj{S - x) dx 

Jo J -00 

+ / dFnid) / aw{x,d\xo)dx / Gx{S,t - d\x + u)dFj{u). 

Jo J —00 J —00 

Hence, after setting 

(19) Bx{S,t\xn) ■.= FR{t)Gw{S,t\xo)+ [ dFnid) [ aw{x,d\xo)Fj{S-x)dx 

Jo J -00 

+ / FR{t-d)dFR{d) I aw{x, d\xo)dx / Gw{S,t ~ d\x + u)dF.j{u) 

Jo J — 00 J — 00 

it is not hard to see that from ^TE\i the following lower bound for the FCT distribution 
function is obtained: 

(20) GxiS,t\xo)>BxiS,t\xo), t > 0. 
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4.1 An improved bound 

Bound (j20p holds in the case of generaUy distributed jumps Ji and renewals Ri. By making 
use of a method already adopted in the proof of Theorem 3.3 of Di Crescenzo and Pellerey 
[lOj . hereafter we see how such bound can be improved by specifying the distributions Fj 
and Fr. Let us assume that 

(i) the random times separating consecutive jumps are exponentially distributed, with 
Ffl(i) = e"^*, t > 0, and A > 0, and that 

(ii) the jumps have fixed amplitude a.s., with Fj(x) = l{x>a} a > 0, which 
corresponds to assumption (|15p with 6 = and 77 = 1. 

In the following we shall stochastically compare {X{t)} with another Wiener process 
with jumps {Xn(t)}t>o, which is driven by the same diffusive component of {X(t)}. We 
formally have 

(21) X„{t)^W{t) + Ynit), t>0, 

where {M^(t)} and {Yn{t)} are independent stochastic processes, {W^(t)} is the same process 
appearing in the right-hand-side of ([T]), {Yn{t)} is a jump process such that Yn{0) = 

and Yn{t) = Y.i=i * > 0' with y„(i) = when N{t) = 0. Moreover, Ji,J2,... are 
real- valued i.i.d. r.v.'s degenerating at na, with a > and n a fixed positive integer, i.e. 
possessing cdf Fj{x) — l{x>ria}- Furthermore, {N{t)}t>o is a counting process independent 
of { Ji, J2, . . .} and characterized by i.i.d. Erlang-distributed renewal times R2, ■ ■ ■ having 
survival function Fj^{t) = e""^* ^^^j ^ > 0, with A > 0. We stress that processes 

{X{t)} and {Xn{t)} share the two parameters A and a. 

Making use of a customary tecnique based on the constructions of "clone" processes (see, 
for instance, Theorem 3.3 of Di Crescenzo et al. '%]) and recalling ^ and pT|) . it is not 
hard to prove that 

(22) Xn{t) <st X{t) for all n > 1, 

where <st denotes the usual stochastic order (see Section l.A of Shaked and Shanthiku- 
mar [2Q]). An immediate consequence of (I22|) is that the FCT's of those processes are 
stochastically ordered too, i.e. 

(23) Tx <st T^^ for all n > 1. 
In other terms, 

(24) Gx{S,t\xo)>G^^{S,t\xa) for all n > 1 and t > 0, 

where G^^ {S, t \ xo) is the FCT cdf of {X,,{t)}. Similarly to ^ the following bound holds 
for all n > 1: 

(25) G^J5,i|xo) >B^J^,t|xo), i>0, 
where 

(26) i3^J5,i|xo) := e-^'Xi^^GH'(5,i|xo) 

+ / Ae-^"^ ^ <ld I aw{x,d\xo)dx 

Jo (n-iy. Jg_^^ 

^ Jo f- {n~l)r 
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/S — na 
aw{x, "i? I xo) Gw{S, t — d\x + no) da;, 
-oo 

with Gw and aw defined in Q and (fT2|) . respectively. In conclusion, by p4|) and ((25l) we 
obtain 

(27) Gjf(S',i I Xo) > sup % (S',i I Xo), t > 0. 

T!>1 

Since (|20|) provides a bound of type Gx{S, t \ xq) > B^^ {S, t \ xo), Eq. ([27]) yields a better 
bound. Furthermore, the right-hand-side of (|27p is not necessarily a distribution function, 
so that we can improve the bound as follows: 

(28) Gjf (S',i I Xo) > G£(S',t I Xo) = max sup Sj^ (5, T I Xo), t>0. 

0<T<t „>1 

Figure 3 shows estimates of FCT distribution function Gx{S,t \ Xq) obtained from sim- 
ulations perfomed by means of the procedure given in Di Crescenzo et al. [7] , together with 
the respective lower bound ([28|l . We stress that the lower bounds given in Figure 3 shows 
a case in which, for all r > 

(29) sup Bj^ [S, r I xo) = B^ (S, r | xo). 

n>l 

Note that ([M)) holds because relation S — na > is satisfied only if n = 1. Instead, the 
angular points appearing in the lower bound on the right of Figure 4 show a case in which 

sup Sjj {S,t\xo) ^ B^ (S,t\xo), 

n>l 

for different values of fc as r varies. 
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Figure 4: On the left: Estimated FCT density gx{S, t \ xq) and lower bound ge{S, t \ xq) for 
a = 2, 6 = 0, xo = 0, S" = 6, A = 0.2, 77 = 1,^^ = 1.5 and cr^ = 0.2. On the right: The 
corresponding estimated distribution function Gx {S, t \ xq) and lower bound Ge{S, t \ xq). 



5 Remarks on computational aspects 

Let us now point out some computational features related to the evaluation of bound 
The last integral in ([^S]) . i.e. 

S — na 

aw{x, 1} I a;o) Gw{S, t — Six + na) dx, 
can be expressed in terms of integrals of the form 
(30) e-''("=-'')' $(ca; + d) da;, 







with a > 0, with 6, c, d e M and $ given in (fTUl) . Unfortunately no closed forms are available 
for (j30p . so that we have been forced to numerically evaluate it by splitting the integration 
domain as follows: 

1=1 e-'^f^^'')' $(ca; + d) dx + / e"'^'^-'*)' $(ca: + d) dx. 

Jo Jb 

Then, transforming the second integral we obtain: 

(31) I = j^e-"^''-^^^ ^{cx + d)Ax + ^^j^ x^''^e-'' ^(^^j^ + ch + d^dx. 

A Gauss-Legendre quadrature rule has been applied to the first integral by choosing 16 
abscissas, after an ad hoc transformation of the integration domain in [—1,1], whereas 
a generalized Gauss-Laguerre quadrature rule has been applied to the second integral by 
choosing 16 abscissas. Note that the latter quadrature formula is exact for weighting func- 
tions as w{x) = x^e~^, with /? S M, which includes the case of the second integral in (|3ip . 

In order to evaluate the goodness of lower bound pT]) we shall now adopt the following 
"measure of closeness" : 

n 

(32) £ = h^[gx{S,jh\xo) ~ ge{S,jh\xo)], 
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1 


0.8 


0.5 


0.2 


0.0 


1.0 


3.9358e-02 


1.1393e-01 


2.4054e-01 


3.6791e-01 


4.5211e-01 


1.5 


2.1059e-02 


6.5103e-02 


1.4572e-01 


2.2688e-01 


2.8409e-01 


2.0 


1.7266e-02 


4.4515e-02 


9.6393e-02 


1.5005e-01 


1.9656e-01 


2.5 


1.7283e-02 


3.4975e-02 


7.0264e-02 


1.1009C-01 


1.4486C-01 



Tabic 1: Values of £ for various choices of /i and rj. 



where h is a, discretization step and n is the minimum integer such that 

n 

hJ29x{S,jh\xo)>0M. 

Differently stated, (15^ measures how the bound gi is close to gx by the multiples of h up 
to the 0.99-percentile. For h = 0.01, Table 1 shows values of £, where is an estimate of 
gx, whereas gi is the lower bound given in (|17p. It is evident that £ increases if fi decreases 
or if T] decreases. 

We finally point out that throughout the paper the estimated functions gx and Gx 
have been obtained by use of 10^ simulated FCT's. The estimated pdf gx has been built 
by adopting an Epanechnikov kernel estimator, while Gx derives by use of a Kaplan-Meier 
estimator. 
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